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The dynamics of charge transfer between molecular species in the condensed 
phase remains a difficult and compelling problem from a theoretical perspective. 
The difficulty stems from the fact that one is faced with the situation that for 
most of the time there is a clear separation in the time-scales between the motion 
of the electronic degrees of freedom and the nuclear or molecular degrees of free- 
dom. Hence, most of the time, nuclear reorganization dynamics following some 
dramatic change in the electronic structure occurs on the Born-Oppcnheimer po- 
tential energy curves parameterized by the nuclear coordinates. However, this 
description breaks down whenever the electronic wave function changes rapidly 
in the direction of the nuclear motion leading to transitions between the dis- 
crete electronic states. This coupling becomes even more dramatic when there 
is an actual crossing of the potential curves at certain nuclear geometries. Even 
more so, the number of nuclear degrees of freedom, while finite, typically far out 
number the number of relevant electronic states for a given process. 

For the case of a infinite thermostat, the Pauli master equation and Redfield 
equation have long been applied to the study of quantum relaxation. While 
originally derived from more or less heuristic arguments, such equations can 
be obtain from formally exact master equations using projection operator tech- 
niques and a series of well defined approximations. In the limit of a continuum 
of modes, the resulting transition probabilities become identical to what one 
obtains using Fermi's golden rule. 

In this paper, we review our recent work in deriving a formally exact and 
time-local approach for incorporating non-Markovian dynamics in to the quan- 
tum master equation for state to state electronic transitions. [l| Here, we discuss 
some of the details of this work with particular attention to its theoretical devel- 
opment. As an application, we examine the kinetics of charge transfer between 
two co-facially stacked conjugated polymer chains using a model recently 

developed by our group. |5j Here we compare the inclusion of both fast and slow 
phonon modes affect the transfer of a hole from one chain to the next following 
photo-excitation. 



The Hamiltonian describing a wide range of photphysical problems can be 
cast in the form: 

H = ^e„|n)(n| + ^g nmi \n)(m\(a\ + a,) + ^c^aja;. (1) 

n nmi i 

Here |n)'s denote electronic states with vertical energies e n , a] and are the 
creation and annihilation operators for the normal mode i with frequency w,, 
and g n mi are the coupling parameters of the electron-phonon interaction which 
we take to be linear in the phonon normal mode displacement coordinate. 

We can separate H into a part that is diagonal with respect to the electronic 
degrees of freedom, 

H = ^e n \n)(n\ + ^2 g nni \n) (n\(a\ + a,-) + ^W;ata i; (2) 

n ni i 

and an off-diagonal part V 

V = ^2 9nmi\n){m\(a\ + en), (3) 

nmi 

where the prime at the summation sign indicates that the terms with n = m are 
excluded. This separation is useful for the following two reasons. First, in many 
systems only off-diagonal coefficients g n mi are small compared to g n m- Hence, V 
can be treated as a perturbation. Second, for many cases of interest, the initial 
density matrix commutes with Hq. In this case, the separation gives simpler 
forms of the master equations. 

For further analysis it is convenient to perform a polaron transform using 

JJ = e 1 /x lv * — > \n)(n\e ^* u * v 1 ; (4) 

n 

in which our transformed Hamiltonian becomes 

H Q = U-'HoU = ^n\n)(n\ + J2^i4 a h (5) 

n i 

where the renormalized electronic energies are 

e„ = e n -y2—- (6) 

% 

Applying the same unitary transformation to V gives 

V = J2\n)(m\M nmi ., (7) 

nmi 
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where the system-bath operators are 

AT ft. 2 9nm \ £ ^ nn} -j, mmj ) ( t } 

M n ™ = # njm I a! + a 4 — le J ^ J . (8) 

At this point it is useful to connect the various terms in our Hamiltonian with 
specific physical parameters. The terms involving (g nn j — gmmj)/^j can be 
related to the reorganization energy 

eL - E (g "" J ~Jr mi? = E 

3 3 3 

where Sj is the Huang-Rhys factor for mode j which is related to the Franck- 
Condon factor describing the overlap between the Vj — 1 vibronic state in one 
electronic state with the Vj = vibronic state in the other. Likewise, the energy 
difference between the renormalized energies is related to the driving force of 
the transition, 

In the transformed picture the electronic transitions from state \n) to |m) are 
accompanied not only by the creation or annihilation of a single phonon of mode 
i but also by the displacements of all the normal modes. This is quite different 
from the spin-boson model which does not have coordinate dependent coupling 
between the electronic states. 

For a properly chosen projection operator V and the initial total density matrix 
that satisfies p(0) = Vp(0), Vp(t) can be shown to satisfy at least two different 
master equations: the Nakajima-Zwanzig (NZ) equation |a, S 13 an d the time- 
convolutionless (CL) master equation 



dVpjt) 

dt 
dV P {t) 

dt 



f drJC NZ (t^r)Vp{r), (9) 
Jo 

[ dTjC CL (T)V P {t). (10) 

Jo 



The explicit expressions for superoperators IC nz (t) and K, cl (t) can be found 
in Ref. [l0| . Since Eq. ^| is less well known than Eq. [5] we will show here how 
Eq. ^|can be derived from Eq. [5] Applying the Laplace transformation to both 
sides of Eq. we obtain 

sVp(s) - Vp(0) = -K,{s)Vp{s), (11) 

where p(s) and /C(s) are the Laplace transforms of p(t) and IC(t), respectively. 
Solving the last equation for p(s) and applying the inverse Laplace transforma- 
tion we have the formal solution of Eq. 

Vp(t) = V(t)Vp(0), (12) 
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where 



K + ZOO 



v v = hL^TTm) (13) 

and k is an arbitrary positive constant chosen so that the contour of integration 
lies to the right of all singularities of the integrand. Differentiating Eq. 1 121 with 
respect to time we obtain the time local master equation 

d ^=nt)v P {t) (14) 

with 

Ht) = d -^v-Ht). (is) 

Eq. ^]can be written in the form of Eq. 1 101 with 

IC CL (r) = -^1. (16) 

Even though Eqs. l9l and 117)1 are formally exact, it is impossible to determine 
K. nz (t) and IC cl (t) for most realistic systems. Since we have assumed the 
coupling to be weak, we ond order in the coupling constants 

K$ Z (j) = IC% l (t) = PC v e^ CoT C v r. (17) 

Here £o and Cv ar e the Liouville superoperators corresponding to Hq and V 
whose action on some density matrix p is given by 

£ oP = Hop - P H n , £ v p = Vp-pV. (18) 

In Ref. ^3 j Breuer and Petruccione show that to second order in the coupling 
constants, the convolutionless experession (Eq. 110(1 gives a better approximation 
to the exact solution than the Nakajima-Zwanzig equation (Eq. EJ). It also has 
an additional mathematical convenience of being local in time. Therefore, in the 
following analysis we will use the convolutionless approach. 

Taking our initial density matrix in the transformed representation as 

p(0) = \n)(n\ 1 . (19) 

we will use the projection operator that acts on the total density matrix in the 
following way 

Vp = J2\n)(n\p^Tr(\n)(n\p), (20) 
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where the trace is taken over both electronic and oscillator degrees of freedom 
and 

P°ea = ( 21 ) 

q Tr(e-^i^ a i a *) 

Note that 

Tr(|n)(n|p) = P n , (22) 

where P n is the probability to hnd the electronic subsystem in state \n). Using 
Eq. the definition of £■§ (t) in Eq. El and the definition of V p given by 
Eq. |20| the following explicit convolutionless equation is obtained 

dP 

m rn 

The time dependent rates W nm (t) are given by 

W mn {t) = 25fte f dTV{M m! M ranj (r)) e - !; ( e "»-^» T , (24) 
Jo ij 

where 

(M nmi M m nj(T)) = Tr {M nml M mnj ( T ) P l s q ) . (25) 

Due to the explicit form of operators M nm i (Eq. the calculation of the corre- 
lation functions in Eq. 1251 can be reduced to the averaging of the displacement 
operators over the equilibrium ensemble (Eq. I21|) . After straightforward, but 
lengthy, calculations we obtain the principal result of our work: 

{■^■nmi-^mnj (^~)) — 9nmi9mnj 

X (y^nmjijlj ~t~ l)c J ^nmjTlj^ J ~t~ ^nmjj 

+ <%(n, + l)e^ r + %n ie -^ T ) 9 nro (r)/ nm (r). (26) 



Here 



(dnni Qrarai) 



(27) 



q (,9nni + Qrami) fORA 

«i = -5-^ . (29) 

?«m(r) = e' E ^ A "- sin ^ T , (30) 

/„ m (r) = e - 2 £ ] ^+D A L J (i-™s S r) > ( 31 ) 
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FIG. 1: (Left) Molecular geometry of TFB:F8BT interface. (Right) Relative energies 
of highest occupied and lowest unoccupied orbitals of TFB and F8BT. 

In the case when all the diagonal electron/phonon terms vanish, g nn i = 0, the 
correlation functions in Eq. |2EI reduce to those obtained within the golden rule 
approach [Tll[T^ |. It is clear from Eq. [25] that for systems with g n ni/^i ^> 1 the 
golden rule approach is not applicable. Unfortunately, the complicated form of 
the correlation functions precludes much further analysis at this point and we 
need to turn to numerical calculations. 

As an application of our approach, we consider the electronic relaxation in 
the conjugated polymer heterojunctions previously investigated in Ref. [Hi , IhH 
by a different approach. Here we consider the decay of an excitonic state into 
a charge-separated state at the interface between two semiconducting polymer 
phases (TFB) and (F8BT). The structure of the polymers at the interface is 
shown in Fig.^ Such materials have been extensively studied for their potential 
in organic light-emitting diodes and organic photovoltaics 0, 0, 0|. At the 
phase boundary, the material forms a type-II semiconductor heterojunction with 
the off-set between the valence bands of the two materials being only slightly 
more than the binding energy of an exciton placed on either the TFB or F8BT 
polymer. As a result, an exciton on the F8BT side will dissociate to form a 
charge-separated (exciplex) state at the interface. @>0>ElEiEl 

F8BT* : TFB — > F8BT~ : TFB+. 

Ordinarily, such type II systems are best suited for photovoltaic rather than LED 
applications. However, LEDs fabricated from phase-segregated 50:50 blends 
of TFB:F8BT give remarkably efficient electroluminescence efficiency due to 
secondary exciton formation in the back-reaction 

F8BT~ : TFB+ — > F8BT* : TFB, 

as thermal equilbrium between the excitonic and charge-transfer states is es- 
tablished. This is evidenced by long-time emission, blue-shifted relative to the 
emission from the exciplex, accounting for nearly 90% of the integrated photo- 
emission. 
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FIG. 2: (Left) The population of the excited state as a function of time for the 
TFB:F8BT heterojunction for three different temperatures obtained by solving the 
time-convolutionless master equation (black) or through the Marcus-type approxima- 
tion (gray) at T = 230K, 298K, and 340K. (Right) Transition rates Wi 2 (t) and W 2 i(t) 
for the TFB:F8BT heterojunction involving 24 normal modes at 298 K. 

We consider only the two lowest electronic levels corresponding to 

\XT) = F8BT* : TFB & \CT) = F8BT~ : TFB+. 

We take the vertical energies for these two states as €qt = 2.191eV and 
exr = 2.294 eV. As was shown in Rcf. [H E |3 E3, in poly-fluorene based 
systems, there are essentially two groups of phonon modes that are coupled 
strongly to the electronic degrees of freedom as evidenced by their presence as 
vibronic features in the vibronic emission spectra, namely: low frequency tor- 
sional modes with frequencies between 90 and 100 cm -1 and higher frequency 
C=C stretching modes with frequencies between 1500 and 1600 cm -1 [l]|. The 
vibronic couplings within the model were determined by comparison between 
the Franck-Condon peaks of the predicted and observed spectra of the system. 

As before, we consider the initial state as being prepared in the \XT) state 
corresponding to photoexcitation of the F8BT polymer. The results for the 
higher electronic level population obtained by numerically solving Eq. I|23|l are 
shown on Fig. The time dependence of coefficients Wi2(i) and Wai(*) is 
shown on Fig. [2 As in the previous example the electronic relaxation does not 
follow a simple exponential pattern. The main difference between this model and 
the previous example is that the electronic relaxation time is of the same order 
as the recurrence time for coefficients W\2{t) and Wziif)- As can be seen from 
Fig. |21 the time dependence of these coefficients has the form of approximately 
constant regions abruptly changed at regular recurrence intervals. This type of 
time dependence makes the electronic relaxation look like a series of exponential 
relaxations with changing rates. We do note, however, that transition rates 
W\2{t) and W2\(t) can become negative as the relaxation proceeds towards the 
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FIG. 3: Exciton population keeping only the high frequency (solid), low frequency, or 
using both high and low frequency phonon modes. 

equilibrium population. As with any approximate method, this can lead to a 
lack of positivity in the populations. 

Note also, that because the relaxation does not obey a simple exponential 
rate law, the initial decay slightly "overshoots" the final equilibrium popula- 
tion. This is most evident in the highest temperature case considered here (T = 
340K). Since the \XT) is also the emissive state, photo-emission (not included 
herein) depletes this population on a nanosecond time scale (radiative lifetime) . 
Population back-transfer from the \CT) to maintain a thermal equilibrium pop- 
ulation then leads to the continuous replenishment of the emissive species such 
that nearly all of the CT population contributes to the formation of secondary 
(regenerated) excitonic states 0, U 0, El ■ 

It is of interest to compare the relaxation dynamics in TFB:F8BT heterojunc- 
tion obtained through the application of the time-convolutionless master equa- 
tion with other approaches. As an example we will consider the Marcus-type 
app roach which is widely used to study electron transfer in chemical systems 
[Hill El El- Based upon the model parameters, the driving force and re- 
organization energy are AE 2 i — 0.073 eV and E\ — 0.092 eV, which places 
the system very close to the barrier-less regime as reported by Morteani, et 
al 0. Because of the coordinate dependency of the coupling, some care must be 
taken in deriving the Marcus rates and we give details of this in the appendix. 
The results for the relaxation dynamics using these rates applied to TFB:F8BT 
hetcrojunction for three different temperatures are shown in Fig. [2] 

It can be seen that apart from a more complicated time dependence the master 
equation approach gives faster relaxation when compared to the Marcus-type 
picture. The discrepancy between the two approaches can be explained by the 
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fact that the Marcus approximation is assumed to be valid when kT 3> hu)i 
for all the normal modes. In the case of the TFB:F8BT heterojunction this 
condition is not satisfied for the higher frequency modes. 

Finally, we consider the contributions of each type of intramolecular motions 
to the charge-transfer process. Recall, our model included two non-overlapping 
bands of intramolecular phonon modes: a low frequency band corresponding to 
the torsional motions of the chain and a high frequency band corresponding to 
C=C bond stretching motions. In Fig. [31 we show the contribution of each band 
to the overall charge-transfer process. If we include only one band or the other, 
very little charge transfer occurs while if we include both bands, transfer occurs 
within the first 0.5ps following excitation. This emphasizes the fact that while 
one set of modes may be bringing the system into and out of regions of strong 
coupling between the XT and CT states (likely the high frequency modes), the 
other set of modes are required to help dissipate the energy associated with 
making the transition. 

In conclusion, we present here a brief overview of a new approach we have de- 
veloped for incorporating non-Markovian dynamics into the calculation of state- 
to-state transition probabilities for electronic systems. The time evolution of 
the electronic populations is, in general, more complex then the one obtained 
with the Pauli master equation and depends on the explicit form of the time 
dependent coefficients. The time-dependence of the rate coefficients introduces 
non-Markovian effects due to the vibrational motion into the electronic pop- 
ulation transfers. What is currently lacking, however, is the inclusion of the 
electronic coherences and the important process of decoherence. Such contri- 
butions are quite important and we are currently extending our approach in 
that direction. While more complete and fully quantum approaches (such as 
MCTDH) permit highly detailed analysis of the dynamics of systems such as 
these, the computational cost associated with these approaches remains quite 
high. 
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